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The continuous time random walk (CTRW) model exhibits a non-ergodic phase when the average 
waiting time diverges. Using an analytical approach for the non-biased and the uniformly biased 
CTRWs, and numerical simulations for the CTRW in a potential field, we obtain the non-ergodic 
properties of the random walk which show strong deviations from Boltzmann-Gibbs theory. We de- 
rive the distribution function of occupation times in a bounded region of space which, in the ergodic 
phase recovers the Boltzmann-Gibbs theory, while in the non-ergodic phase yields a generalized 
non-ergodic statistical law. 
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The ergodic hypothesis is a cornerstone of statistical 
mechanics. It states that ensemble averages and time av- 
erages are equal in the limit of infinite measurement time. 
Starting with the work of Bouchaud, there is growing in- 
terest in stochastic ergodicity breaking (SEB) which finds 
applications in a wide range of physical systems: phe- 
nomenological models of glasses Q , laser cooling 0, , 
blinking quantum dots [dja] , and models of atomic trans- 
port in optical lattices [y]- SEB is found for systems, 
whose dynamics is characterized by power law sojourn 
times, with infinite average waiting times. In such sys- 
tems the microscopical time scale diverges, for example 
the average trapping time of an atom in the theory of 
laser cooling |2, l| ■ The relation between SEB and di- 
verging sojourn times can be briefly explained, by noting 
that one condition to obtain ergodicity is that the mea- 
surement time is long, compared with the characteristic 
time scale of the problem. However this condition is never 
fulfilled if the microscopical time scale, i.e. the average 
trapping time, is infinite. It is important to note that 
the concept of trapping time probability density function 
(PDF) tp(t), with diverging first moment, is wide spread 
and found in many fields of physics 0, 0, S @ ■ It was 
introduced into physics by Scher and Montroll in the con- 
text of continuous time random walk (CTRW) 10]. This 
well known model 0, 01 exhibits anomalous sub-diffusion 
(r 2 ) ~ t a with a < 1, and aging behaviors which are 
related to SEB. 

Clearly if the CTRW is non-ergodic, Boltzmann-Gibbs 
statistics is not valid, in a way defined precisely later. 
The goal of this letter is to obtain a generalization of 
Boltzmann-Gibbs statistical mechanics, for CTRW mod- 
els. Besides its theoretical importance this goal is timely 
due to recent observations on the single particle level of 
CTRW type of dynamics 0,0], for example anomalous 
diffusion of a single magnetic bead in a p olymer network 
with a well defined temperature T In single par- 

ticle experiments, the many particle averaging, i.e. the 
problem of ensemble averaging, is removed |14j. Hence 
a fundamental question is whether time averages of sin- 
gle particle trajectories yield information identical to en- 
semble averages. The large number of applications of the 



CTRW model, and related models like the trap model 
and the comb model, make us believe that constructing 
a general non-ergodic theory for such systems is worthy. 

Before introducing the model, recall that the basic tool 
in statistical mechanics is Boltzmann's probability Pb(x) 
of finding a system in a state with energy E(x), using the 
canonical ensemble 



Pb(x) 



exp 



f g(x) 



(1) 



where T is the temperature and Z = exp[— E(x)/T] 
is the normalizing partition function. Here for simplicity 
we assume a discrete energy spectrum. To obtain the 
ensemble average of a physical observable, for example 
the energy of the system, we use (E) = J2 x E(x)Pb(x), 
and similarly for other physical observables like entropy, 
free energy etc. When measurement of a single system 
is made, a time average is recorded. Consider a system 
randomly changing between its energy states {E(x)}. Let 
t x be the total time spent by the system in energy state 
E(x), within the total observation interval (0, t). We 
define the occupation fraction 

p x = — , and the time average energy is E = E(x)p x 

X 

According to statistical mechanics, once the ergodic hy- 
pothesis is satisfied, and within the canonical formalism 
p x = Pb(x) and then E = (E), and similarly for other 
physical observables. Thus for ergodic systems the frac- 
tion of occupation time is non-random, in the thermo- 
dynamic limit of long measurement time 15]. For non- 
ergodic systems the occupation fraction p x is a random 
variable, even in the long time limit. Thus, an important 
goal of the theory of SEB is to calculate the distribu- 
tion function of the random variable p x . We will show 
below, in the context of CTRW models, that a rather gen- 
eral distribution function, describes statistical properties 
of the fraction of occupation times in the non-ergodic 
phase. And that this distribution function is related to 
the partition function of the problem. 



2 



0.4 



o.i 



L J h 




FIG. 1: Boltzmann's equilibrium for an ensemble of CTRW 
particles, in an harmonic potential field, and fixed tempera- 
ture T = 3. In simulations (cross) we use a = 0.3, 0.5, 0.8, 
the results being indistinguishable. The figure illustrates that 
for an ensemble of particles standard equilibrium is obtained. 
Ergodicity breaking is found only when long time averages of 
single particle trajectories are analyzed. The scaled poten- 
tial (dot dash curve) is the harmonic potential field, and the 
theoretical curve is Boltzmann's equilibrium distribution. 



We consider a one dimensional CTRW on a lattice. 
The lattice points are labeled with index x and x — 
—L, —L + l, 0, ...L, hence the system size is 2L+1. On 
each lattice point we define a probability < Qr(x) < 1 
for jumping right, and a probability for jumping left 
Ql(x) = 1 - Q R (x). Let be the PDF of waiting 
times at the sites. The particle starts at site x = 0, it 
will wait there for a period t\ determined from i/>(i), it 
will then jump with probability Ql(0) to the left, and 
with probability Qr(0) to the right. After the jump, say 
to lattice point 1, the particle will pause for a period i2, 
whose statistical properties is determined by tp{t). It will 
then jump either back to point x = or to x = 2, ac- 
cording to the probability law Qr{1). Then the process 
is renewed. We assume reflecting boundary conditions, 
namely Ql{L) = Qr{—L) = 1. We consider the generic 
case 8, 9], where 



\T(-a)\ 



(2) 



when t — > oo and < a < 1, A > 0. Specific values of a 
for a wide r ang e of physical systems and models are given 
in 0, H 0, llOl IrH flaj . In this case the average waiting 
time is infinite. 

We obtain now an important limiting distribution 
function, which will yield the statistical properties of the 
non-ergodic CTRW. We consider a specific lattice point 
x. We introduce a state function 9 x (t) which is equal 
1 when the particle is on x, other wise it is zero. Thus 
9 x (t) jumps between the value 1 (state +) and zero (state 
— ) and vice versa. The PDF of times when the particle 
occupies state + [— ] is denoted with ip+(t) [ip-(t)] re- 
spectively. In our model ip+(t) — ip(t). To obtain ip~(t) 



note that after the particle leaves lattice point x it is ei- 
ther on x + 1 or on x — 1. Let tji (t^) be the random time 
it takes the particle starting on x + I (x — 1) to return 
to x, and /r^r) [/l(^l)] the corresponding PDF of the 
first passage time respectively. Then the PDF of times 
in state — is 



V>- (t) = Q R (x)f R (t) + Q L {x)f L (t) . 



(3) 



We will later find explicit expressions for ip_(t), however 
for the time being let us assume that in the limit of long 
t it behaves like 



<Mt) 



|r(-a)| ' 



(4) 



where A x > will depend on model parameters. Let 
t x be the total time spent on point x, within the time 
period (0, t). Then the occupation fraction is p x — t x /t — 
Jq 0x(t')dt' /t. A calculation, whose details will be given 
elsewhere, shows that the PDF of the occupation fraction 
in the limit of infinite measurement time is 



f(Px) = $aCR.x,P x ), 

where 1Z X — A/A x and 5 a {TL x ,p x ) = 



(5) 



sin ira 



n x r x -"{i-p x y 



7T TZ X 2 (1- p x y a +p 2 x a + 2TZ X {1 -v x Tv a x cos to 

(6) 

This equation indicates SEB since the occupation frac- 
tion remains a random variable, even in the limit of long 
measurement times. The PDF Eq. (@ is normalized ac- 
cording to f Q 5 a (lZ x ,p x )dp x = 1. When TZ X = 1, Eq. 
© was obtained by Lamperti 0] in the context of the 
mathematical theory of occupation times 0|j an d see 
5] for a physical application in the context of blinking 
quantum dots. In particular when 1Z X = 1 and a = 1/2 
we find the arcsine distribution. More generally the am- 
plitude ratio 1Z X determines the degree of symmetry in 
the problem as we will demonstrate later. 

We use a general physical argument to obtain 1Z X . As- 
sume that the random walker is in contact with a thermal 
heat bath, with temperature T, and interacting with an 
external potential field E(x) . For an ensemble of particles 
Boltzmann Gibbs statistics must hold. In particular the 
probability that a single member of an ensemble of parti- 
cles, will occupy the lattice point x, is given by Pb{x) Eq. 
(Q. Eq. |jfj| shows that p x is a random variable, how- 
ever when we ensemble average the occupation fraction 
we must obtain Boltzmann-equilibrium statistics 



(Px) = / Pxf (Px) dPx = Pb{x). 
Jo 

On the other hand Eq. © yields 
,_ , Tlx 

\Px 



(7) 



(8) 
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FIG. 2: The PDF of occupation fraction p x — t x /t where t x 
is the occupation time on lattice point x — 0. The random 
walk is in an harmonic potential field, the point x — being 
the minimum of energy. For an ergodic process satisfying 
detailed balance, the PDF f(p x ) would be narrowly centered 
around the value predicted by Boltzmann which is given by 
the arrow. In a given numerical experiment, it is unlikely 
to obtain the value of p x predicted using Boltzmann-Gibbs 
ergodic theory. The solid curve is the analytical formula Eq. 
<tTlt with a = 0.3 and T = 3 



Hence we find 



Pb{x) 
1-P b (x) 



exp 



(9) 



where Z' = J2 y ex P(~ E(y)/T) and the sum is over all 
energy states excluding state x, Eqs. (0 describe 
the relation between the non-ergodic dynamics and the 
partition function of the problem. Thus we established 
an explicit relation between SEB and the basic tool of 
equilibrium statistical mechanics. The remaining goal of 
this Letter is to prove our physical picture, based on the 
CTRW model. 

We now consider the unbiased one dimensional CTRW 
on a lattice x = —L, with Ql(%) — Qr(x) = 
1/2. Such a process describes sub-diffusive motion. We 
obtained the long time behaviors of PDFs /«(<;) and /z,(i) 
(details are left for longer publication) , and then obtained 
the PDF of the fraction of occupation time p x on a lattice 
point x, excluding the boundary points, 



lim f(p x )=6 a ((2L-l)-\p Q 



(10) 



Eq. HI U| l does not depend on the position x of the obser- 
vation point, reflecting the symmetry of the problem. 

The biased CTRW yields anomalous diffusion with a 
drift. For this case the probability of jumping left is 
Ql{x) = q, and right is 1 — q where q ^ 1/2. Unlike the 
unbiased case now clearly different locations along the 
lattice, have different distributions of the fraction of oc- 
cupation time. A detailed calculation of the first passage 



FIG. 3: Same as Fig. J2J however now a = 0.8. Instead of 
the U shape found in Fig. we find a distorted W shape 
of the PDF. A peak close to Boltzmann's value for p x , i.e the 
arrow on Pb(x) is an indication that as a is increased the 
ergodic phase is approached. 



times, shows that in the long time limit, 



(11) 



with 1Z X 



2q-l 



1-9 



-(i-«r 



1-9 

q 



- 1 



(12) 

Note that Eq. pi(l is not sensitive to the short time 
behavior of the waiting time distribution. When q = 1/2 
we recover Eq. lfTU)l. 

The biased CTRW is used to model anomalous diffu- 
sion under the influence of a constant external driving 
force e.g. t 10]. If the physical process is close to 
thermal equilibrium, the condition of detailed balance is 
imposed on the dynamics. This standard condition im- 
plies that for an ensemble of particles Boltzmann's equi- 
librium is obtained. The potential energy at each point 
x, excluding the reflecting boundaries, due to the inter- 
action with the external driving force is E(x) = —Tax 
and a is the lattice spacing. The condition of detailed 
balance then reads 



1 



1 + exp ' 



(13) 



Using Eqs. (|12I13|I we can rewrite the solution in an 
elegant form 



f(Px) = S a 



Pb{x) _ 
\-P B {xY K 



(14) 



where Pb (x) is the canonical Boltzmann factor Eq. Q). 
When the external force is zero we have Pb (x) = l/Z and 
Z = 2L. Eq. (|14J) proves that our physical arguments 
leading to Eq. © are valid, at-least for the uniformly 
biased and unbiased random walks. More generally, for 
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CTRWs far from thermal equilibrium, Eq. 1 (14(1 is still 
valid however one must replace Pb(x) with the corre- 
sponding equilibrium probability. 

Eq. (|14|) shows that the fluctuations p x are vast, in 
particular in any single measurement we are not likely 
to measure the averaged value Pb (x) (see details below) . 
In the limit a — > 1 our theory reduces to the standard 
canonical theory, since then / (p x ) = 5\p x —Pb(x)]. Using 
Eq. HI 4(1 we can explain the meaning of the asymmetry 
parameter 1Z X in Eq. Q. The numerator of 1Z X is Pb(x), 
namely the probability (in ensemble sense) of finding the 
system in state + (i.e. 9 x (t) = 1) while the denominator 
1 — Pb (x) is the probability of finding the particle in state 
— (i.e. 6 x (t) — 0). This interpretation suggests to us that 
Eq. 1 (14(1 has a more general validity, beyond the biased 
and non-biased CTRW. 

More generally we define an energy profile for the sys- 
tem {E-l,E_l + i, E x , ■ ■ •}. The well known detailed 
balance, relates between density of particles at equilib- 
rium N eq (x) at points x and say x — 1 according to 

N cq {x)Q L {x) = N cq {x - 1) [1 - Ql{x - 1)] (15) 

where N cq (x)/N cq (x — 1) = exp f S^-zEz ^ anc [ hence 
we get the constrain on the transition probabilities 
Ql(x)/[1 - Q L (x - 1)] = exp (- E -%T E ' y For a gen- 
eral energy field we postulate that, if the energy profile 
yields a Boltzmann-Gibbs ergodic behavior for a waiting 
time distribution with finite moments (e.g., exponential 
waiting times), then for the same energy profile and a 
long tailed waiting time PDF ip(t) given in Eq. (J2J, our 
central Eq. ((14(1 is still valid. Now Pb(x) depends of 
course on the specific energy profile under investigation. 

We check numerically the generality of Eq. 1 (14(1 , using 
the example of a random walk in an harmonic poten- 
tial. The problem of anomalous diffusion in harmonic 
field was considered in the context of fractional Fokker- 
Planck equations [lij and in single protein experiments 
[l2|. As a byproduct, our work shows that fractional 
Fokker-Planck equations can be used to describe den- 
sity of many particles and not time average quantities, in 
this sense the fractional kinetic framework is very differ- 
ent than the standard Fokker-Planck equations. 

The potential field we choose is E(x) = Kx 2 , with 
K = 1, and T = 3. Eq. 1 (15(1 and the symmetry con- 
dition Ql(0) = 1/2 define the set of transition proba- 
bilities {Ql(x)} for the problem. In Fig. ^ we check 
that our simulations yield Boltzmann equilibrium in the 
Harmonic field for an ensemble of particles. We then 
consider one trajectory at a time. We obtain from the 
simulations, the total time t x spent by a particle on the 
lattice point x — 0, namely at the minimum of the po- 
tential, and then construct histograms of the occupation 
fraction p x =t x /t. 

We consider the case a = 0.3 in Fig. an d show an 
excellent agreement between our non-ergodic theory Eq. 



((14(1 and numerical simulations. The figure exhibits a 
U shaped PDF. To understand this behavior, note that 
for a << 1 we expect that the particle will get stuck 
on one lattice point during a very long period, which is 
of the order of the measurement time t. This trapping 
point, can be either the point of observation (e.g. x = 
in our simulations) or some other lattice point. In these 
cases we expect to find p x ~ 1 or ^ ~ 0, respectively. 
Hence the PDF of p x has a U shape. This is a strong 
non ergodic behavior, in the sense that we have a very 
small probability for finding occupation fraction close to 
the value predicted based on Boltzmann's ergodic theory 
(the arrow). 

When we increase a we anticipate a "more ergodic" 
behavior, in particular in the limit a — > 1. An ergodic 
behavior means that the PDF of the occupation fraction 
p x is centered on the Boltzmann's probability (i.e. the 
arrows in the Figs). In Fig. 0we set a = 0.8 and ob- 
serve a peak in the PDF of p x centered in the vicinity of 
the ensemble average value. Note however that the PDF 
f(p x ) still attains its maximum on p x = and p x = 1. 

More generally, Eqs. (0 |5J) yield the non-ergodic sta- 
tistical mechanical theory of the CTRW model, in the 
sense that our theory gives the distribution of p x , while 
the ergodic Boltzmann-Gibbs theory states p x = Pb (x) . 
Our arguments leading to Eqs. 10 © are general, hence 
our theory might not be limited to CTRW models. The 
mathematical foundation of the theory is the limit the- 
orem (0 related to the arcsine law. The physical input 
is the anomalous diffusion exponent a. A connection be- 
tween the non-ergodic dynamics and the partition func- 
tion was found, which enables us to find non-trivial SEB 
properties of the underlying random walk, in particular 
the random walk in a potential field. 
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